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Strong coupling expansion is computed for the Einstein equations in vacuum in the 
Arnowitt-Deser-Misner (ADM) formalism. The series is given by the duality principle in 
perturbation theory as presented in [M.Frasca, Phys. Rev. A 58, 3439 (1998)]. An exam- 
ple of application is also given for a two-dimensional model of gravity expressed through 
the Liouville equation showing that the expansion is not trivial and consistent with the 
exact solution, in agreement with the general analysis. Application to the Einstein equa- 
tions in vacuum in the ADM formalism shows that the spacetime near singularities is 
driven by space homogeneous equations. 
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1. Introduction 

Strong coupling expansions in theoretical physics have a longstanding tradition. One 
should consider fluid dynamics where two different regimes of motions were found 
depending on what term of the Navier-Stokes equation is taken as a perturbation ^ . 
Boundary layer problems in perturbation theory firstly appeared in this way ^. In 
the seventies attempts to build a strong coupling expansion for quantum field theory 
were put forward by Bender and colleagues 3''*'^. These latter expansions faced the 
problem of the dependence on a cut-off that was not possible to dispose of in a 
well definite manner. The idea behind was simple and smart: For a self-interacting 
scalar field the perturbation is not the potential but the kinetic term. 

What makes these methods rather interesting is that they exploit an evident 
symmetry in the choice of the perturbation in a differential equation. This choice is 
arbitrary and different perturbation schemes can be devised depending on it. So, one 
may ask what is the relationship between different series. This matter was exploited 
in Ref.6 where a theorem was proved showing that in the strong coupling expansion 
for a quantum system the leading order is obtained by the adiabatic theorem in 
quantum mechanics. This striking result was further analyzed in Ref.7 showing how 
the proper expansions can be obtained in different regimes for some Hamiltonians 
in quantum optics. A duality principle holds that interchanging the terms in the 
perturbation series means to obtain a series with the reversed ordering parameter. 
A successful application of this idea in quantum field theory has been recently given 

8 



1, 2008 21:33 WSPC/INSTRUCTION FILE LANL 
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The leading order is given by the adiabatic theorem in quantum mechanics for 
the Schrodirigcr equation but this result is indeed ubiquitous and we will show it 
as a generic result for a strong coupling expansion for a differential equation. This 
theorem will permit us to extend the above analysis to Einstein equations in vacuum 
in the form given in the Arnowitt-Deser-Misner (ADM) formalism ^. Indeed, a 
strongly perturbed system, understood as a system undergoing a largely increasing 
perturbation, becomes a system with an increasingly slow dynamics, that is, an 
adiabatic system. More importantly for our aims is the fact that these equations 
behave homogeneously with respect to the spatial variable in a regime that appears 
dual, in the sense of perturbation theory ^, to the case of the gravitational waves 
that represent small perturbations on a given spacetime. 

The relevance of this result for the Einstein equations is easily understood being 
open up several problems that presently only a numerical treatment can manage 
with all the shortcomings that this implies lo^^. One of the questions that could 
benefit from a strong coupling analysis is the singularity problem Indeed, it is 
well-known that while theorems exist due to Penrose and Hawking I3,i4,i5 pj-gving 
the existence of singularities in the solutions of the Einstein equations, but very few 
is known about the nature of such solutions and their behavior near the singularity 
itself. Our approach permits to analyze the Einstein equations in a regime of very 
intense gravitational field and then to get an answer to this problem. Indeed, our 
approach permits to conclude that the dynamics near the singularity is described 
by homogeneous solutions of the Einstein equations and some numerical results hint 
in this direction i^-^^^ 

What we are going to show is that a dual expansion exists to the weak field 
one that is commonly used to study gravitational radiation-'^^. The main result we 
obtain in this way is that the leading order equations are just the homogeneous 
part of the Einstein equations. This reduces the spatial dependence to parameters 
in the integration constants and the gauge choice at this order. This in turn means 
that in a regime of strong gravitational field the solutions to be considered are the 
homogeneous ones. We will see that in two dimensions the situation is quite similar 
producing a non-singular solution but with an analogous duality between weak and 
strong field solutions. As a by-product we will obtain a strong coupling expansion for 
the gravity in two dimensions as described by the Liouville equation^*'^^'^". We will 
show the agreement with the exact solution giving full consistency to our approach 
for a well known non-linear case. This will yield a start up to our perturbation 
procedure toward the final proof. 

The paper is so structured. In sec. 2 we present duality in perturbation theory 
and the adiabatic approximation at the leading order. In sec. 3 the approach, in a 
general ADM-like setting, is applied to two dimensional gravity. In sec. 4 the method 
is applied to Einstein equations in vacuum in the ADM formalism. In sec. 5 we give 
some applications showing the 1/G behavior of the perturbation series. Finally, in 
sec. 6 conclusions are given. 
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2. Duality in perturbation theory and adiabatic approximation 

In this section we will show how a duality principle holds in perturbation theory 
showing how to derive a strong coupling expansion with the leading order ruled 
by an adiabatic dynamics in order to study the evolution of a physical system. We 
consider the following perturbation problem 

dtu = L{u) + XV{u) (1) 

being A an arbitrary ordering parameter: As is well known an expansion parameter 
is obtained by the computation of the series itself. The standard approach assume 
the limit A — > and putting 

u = uo + Xui + . . . (2) 

one gets the equations for the series 

dtuo = L{uq) (3) 
dtui = L'{uo)ui + V{uo) 

where a derivative with respect to the ordering parameter is indicated by a prime. 
We recognize here a conventional small perturbation theory as it should be. But 
the ordering parameter is just a conventional matter and so one may ask what 
docs it mean to consider L{u) as a perturbation instead with respect to the same 
parameter. Indeed one formally could write the set of equations 

dtvo = V{vo) (4) 
dtvi = V'{vo)vi + L{vo) 

where we have interchanged L{u) and V{u) and renamed the solution as v. The 
question to be answered is what is the expansion parameter now and what derivative 
the prime means. To answer this question we rescale the time variable as t = Xt 
into eq.(l) obtaining the equation 

XdrU = L{u) + XV{u) (5) 

and let us introduce the small parameter e = j. It easy to see that applying again 

the small pertTirbation theory to the parameter e ^ we get the set of equations 
(4) but now the time is scaled as t/e, that is, at the leading order the development 
parameter of the series will enter into the scale of the time evolution producing a 
proper slowing down ruled by the equation 

edfVo = V{vq) (6) 

that we can recognize as an equation for adiabatic evolution that in the proper limit 
e ^ will give the static solution V{uo) = 0. We never assume this latter solution 
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but rather we will study the evolution of eq.(6). Finally, the proof is complete as 
we have obtained a dual series 

u = vo + jvi + ... (7) 

by simply interchanging the terms for doing perturbation theory. This is a strong 
coupling expansion holding in the limit A ^ OC' dual to the small perturbation 
theory A ^ we started with and having an adiabatic equation at the leading 
order. 

It is interesting to note that, for a partial differential equation, we can be forced 
into a homogeneous equation because, generally, if we require also a scaling into 
space variables we gain no knowledge at all on the evolution of a physical system. 
On the other side, requiring a scaling on the space variables and not on the time 
variable will wash away any evolution of the system. So, on most physical systems 
a strong perturbation means also a homogeneous solution but this is not a general 
rule. As an example one should consider fluid dynamics where two regimes dual 
each other can be found depending if it is the Eulerian or the Navier-Stokes term 
to prevail. In general relativity things stay in a way to get a homogeneous equation 
at the leading order. The reason for this is that products of derivatives or second 
order derivatives in space coordinates are the only elements forming the Einstein 
tensor beside time dependence. 

This completes the analysis initially realized in quantum mechanics in Ref.6, 7 
permitting a direct application to general relativity. 



3. Strong coupling expansion for two-dimensional gravity 

As is well-known a direct computation of the Ricci tensor in two dimensions gives 

zero proving that two-dimensional general relativity is trivial. This is due to confor- 
mal invariance of the theory in two dimensions. Different proposals have appeared 
to have non trivial solutions and one of the most common theories is off'ered by the 
Liouville model I8,i9,20 ^j^g^^ j-^^^g ^g^^ applications also in string theory^^. Our 

aim in this section is to obtain a dual perturbation series to the weak perturbation 
series as given in for this model of two dimensional gravity in order just to give 
a scheme of working of the machinery to be also applied to the full four dimensional 
case. So, kit us write the general line element in two dimensions in the form proper 
to ADM formalism 

ds^ = [-a{x,tf + Pi{x,tf]dt'^ + 2l3i{x,t)dtdx + ^{x,t)dx^ (8) 

where we have set a{x,t) for the lapse function, Pi{x,t) for the shift function and 

7(x, t) the metric factor in one dimension. We can simplify the metric by taking the 
shift factor to be zero, so in the following j3i{x,t) = and this gives 



ds"^ = -a{x, tfdi^ + 7(a;, t)dx^ 



(9) 
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that for the equation i? = A, with R the Ricci scalar and A a constant, yields the 
ADM equations 

da = -2aK (10) 

dtK = --K^ - dla + -^d^ad^j + ^aj 
7 27 4 

with the freedom to fix the gauge fimction a. to compute K, the extrinsic curvature, 
and 7. A natural choice in this case is to put a-^ = 7 to make explicit the conformal 
invariance of the metric but other choices are also possible. This gives back the 
Liouville equation 18,19,20 

□0 + Ae'^ = O (11) 

with the introduction of the field cj) = Iwy and being A the only parameter of the 
theory. Indeed, a small perturbation theory is obtained by the explicit solution^^ 

(j){x,t) = (l)o{3^,t) - 2ln(l + ^ J dxidtiG{x - xi,t - ti)e'''^'''^^'^^ (12) 

being <j>o a solution of the equation O(f>o = and G can be expressed by the functions 
g{z) = 9{z) or g{z) = 0{—z) with light cone coordinates as G(x — y) = \9{x'^ — 
y^)g{x^ — y^). As one can see in this case the start-up solution is a solution of the 
free theory. 

Eqs.(lO) can be immediately used with the perturbation method of sec. 2 but 
this application is not so straightforward as one could think. There is a subtle point 
to take into account as also will happen for the full dimensional case. Indeed, we 
note that the lapse function a and the metric factor 7 are both without physical 
dimensions so the extrinsic curvature K has dimension [i] as the derivative. This 
means that when we rescale time variable, we have also to rescale K properly. As 
the variable of our expansion is A we see that we have to do the following scale 
changing 

t \/Ai (13) 
K VAK 

but we can take into account of the rescaling in the extrinsic curvature in its series 
development as 



(14) 



As said we also assume the gauge to be fixed with the conformal choice = 7. 
With the substitution tt = —Kj^ the ADM equations become 

dtOL = aw (15) 
2 \ a 



February 1, 2008 21:33 WSPC/INSTRUCTION FILE LANL 



6 Marco Frasca 



from which is straightforwardly obtained the LiouviUe equation taking into account 
that a = ei . The LiouviUe equation is exactly solvable by two arbitrary functions 
that depend on light cone coordinates 1^,19, 20 ^^gg follows the following set of 
perturbation equations 



aoTTo 

aoTTi + aiTTo 
ao7r2 + OL\-K\ 



(16) 



drTT2 



-aotti + dx 



1 



ao 



-aoa2 - + 



being r = \/Ai and we have adopted for tt the same scaling of the extrinsic curvature 

K and wc have considered a ^ ao + j^cti + ■ ■ ■■ At the leading order these equations 
give the homogeneous LiouviUe equation as it should be expected 

d^cl^o + = (17) 

that has the solution 



/>o(t) = In 



1 



2C| cosh^ 



(18) 



SO that ao = e~ . Ci and C2 are two integration constants that can depend on 
X but not on A otherwise the perturbation series would lose meaning. When C\ 
and C2 are taken not depending on x eq.(18) is an exact solution to the LiouviUe 
equation as can also be proved taking the general solution to this equation and 

expressing it by light cone coordinates. Then, higher order corrections are not zero 
when e.g. Ci is taken to be dependent on x. Assuming for the sake of simplicity 
Ci = kx, being k a constant, the equation to solve is 

<Py , ,'y + A-^/2 

dr2 



2C|cosh^(M) 







(19) 



f 1 — ^0 o 

being y = e 2 whose solution is trivially given by y = — fc that is 



ai 



*1 fc^ *0 

p 2 = p 2 

2 



and to second order 



a2 = e 



^2 



k io 
-— e 2 . 



(20) 



(21) 
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We sec that one has a series in fc^/A as it should be in agreement with the exact 
solution given by 



a = \il — —e 2 



(22) 



and at the end of computation we can fix the constant C2 as to have the scale factor 
^1 — ^ as a renormalization factor giving again an exact solution 



(j){x,T) = -2 In 



cosh 




(23) 



We see that the ground state limit — > — 00 corresponds to the dispersion relation 
k'^ = A producing a null metric. Besides, the dual weak coupling limit k^ ^ A gives 
a "gravitational wave" -like solution in a complete analogy with the four dimensional 
case as we will see. 

Finally, we have effectively got a series expansion in the parameter dual to 
the small perturbation theory and in agreement with our adiabatic hypothesis in 
sec. 2. The leading order is ruled by the homogeneuos Liouville equation. Besides, 
our approach led us directly to exact solutions of the equation itself. 

Some interesting conclusions can be drawn from the "cosmological" solution 
(18). The metric factor is given by 



7(i) 



1 



2C|cosh^(^^) 



(24) 



that has the property to become a constant in the limit r ^ and this situation is 
quite different from the full dimensional theory where from Kasner solution one can 
see that metric blows up. This conformal factor bc;c;omes zero in the opposite limit 
T 00 but generally our strong coupling expansion should have a proper time scale 
given by the rescaling parameter A. As this parameter becomes larger also the time 
scale we are able to analyze increases. In the present case we are in the favorable 
situation that the perturbation terms contain global informations. Finally let us 
point out that such a kind of cosmological solutions was found in string theory . 



4. Strong coupling expansion for Einstein equations 

Our approach for the case of Einstein equations in vacuum will be on the same lines 
as for the two dimensional case. We would like to point out that the ground to reach 
our aim is derived from numerical general relativity as devised e.g. in Ref.23. This 
means that we consider Einstein equations in the ADM formalism in a coordinate 
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basis explicitly given by the set 

dt7ij - P^da.j = ji,d,0' + -ladjP' - 2aK,, (25) 



for a metric 



+7"P [{di-fjn + dj-fin - dn-fij)dl'yrnp 

-didja + ^7'™(5i7jm + dj^im - dmlij)dia 



ds'^ = (-a^ + j,j(3'P^)dt^ + 2(5'dx,dt + j.jdx^dx^ (26) 



being a the lapse function, the shift vector and 7^ the spatial part of the metric. 
The gauge freedom corresponds to fixing a and /?' while jij is to be considered a 
fundamental variable of the theory. Kij is the extrinsic curvature also a fundamental 
variable of the theory to be solved for. Four constraint equations also hold 

^^^R + K"^ - K,jK'^ ^0 (27) 

being ^^^R the Ricci scalar for 7*-' and Vj the corresponding covariant derivative. 
We do not worry about the constraint equations in the following while these are 
relevant matter in numerical computations. 

On the basis of the duality principle in perturbation theory discussed above, wc 
want to prove that in a strong gravitational field the spatial part of the Einstein 
equations, at the leading order, can be neglected and this behavior is exactly dual to 
the case of the gravitational radiation in weak field expansion, a situation modeled 
in a very simple form in the case of the Liouville equation. 

Our analysis starts from a small perturbation theory that can be derived from 
the ADM equations by fixing the gauge freedom as e.g. 

a = 1 (28) 
/?' =0. 

Then we take 

7«=%-+A7g^++A%g' + ... (29) 

being rjij = Sij the flat spatial metric, but other choices are also possible, and in 
this case K^j^ = 0, and A an ordering parameter to define the expansion The 
ordering parameter A is introduced here just as a computational tool to obtain the 
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perturbation scries. As said in sec. 2 this is a standard approach in perturbation 
theory. The proper parameter emerges naturally after the series is computed if the 
series exists and this is the case as we are going to prove both for A — > 0, the known 
case, and A ^ oo. So, the equations we get are 



at7g^ = -2K\f (30) 

^ ij 2 ^-^ 2 ■'^ ' 2 2 
-2/f«iff)+if(i)if« 



that as expected display at the leading order the dynamics of gravitational waves. 
Different choices of the gauge freedom do not change the physics. The leading order 
is ruled just by the linear part of Einstein equations. In this way we recognize the 
standard weak field expansion that holds for A ^ 0. 

We now compute the dual series to the weak field expansion and we will show 
that the situation is quite similar to the two dimensional case with the notable 
difference that solutions of Einstein equations have singularities. As said, in order 
to get a non trivial expansion we have to consider the following settings 

T=\/Xt (31) 

_ (0) I 1 (1) I J_ (2) , 

1 1 

a = ao + —ai + -to Q!2 + • • • 
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while the shift vector /3* is taken to be zero. The following set of equations is so 
obtained 

a.7f = -2aoKlf (32) 



(0) 



(0) 
mi 



._^"P(0) 
2 ' 



;(a.7j°^ + " a„7if )ap7^:^ + a.7l°^5,7(?) ]} 



So, we succeeded in the derivation of a non trivial set of equations for a strong 
coupling expansion of the Einstein equation. We see that in this case we have to 
take A ^ oo to make the expansion meaningful. This is the main result of the paper 
that permits us to draw the following conclusions. At the leading order the set of 
equations that rules the dynamics is that part of the Einstein equations that holds 
for an homogeneous spacetime. So, unless boundary or initial conditions happen to 
depend on spatial coordinates, only the time variable appear to be significant in a 
strong gravitational field. It is interesting to notice here that in the above equations 
the nonlinear part of Einstein equations has taken the role of the linear part that 
we have seen to rule the weak field perturbation theory and a direct evidence of the 
duality principle in perturbation theory as discussed in sec. 2. The situation is exactly 
mirroring the case of the two dimensional gravity by the Liouville equation with 
the difference that here we have a singular behavior. Indeed, we see that due to the 
peculiar structure of the Einstein equations this should not come out as a surprise. 
The most important conclusion to be drawn from this result is that in a strong 
gravitational field, that is near a singularity, the spacetime is indeed homogeneous 
reaching the conclusion we aimed to. This is the well-known Belinski-Khalatnikov- 
Lifshitz (BKL) conjecture 24,25,26 ijn^gg^i ^ one can complete the proof by assuming 
a = ao = 1 in the identical way as we did for the weak field perturbation theory. 
This just simplifies the equations without changing the conclusions consistently with 
the duality principle in perturbation theory. Finally we note that the computation 
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of higher order corrections can be done near a singularity displaying the corrections 
due to irregularities in the spacetime as one departs from the singularity when 
proper boundary conditions are given. 

Two main points should be addressed to complete our proof. Indeed, one should 
understand how general are the above equations with respect to a simple rephras- 
ing of the BKL solutions. But eqs.(32) are able to express more general physical 
situations as also we will see in the next section. A typical situation is given by a 
strongly perturbed black hole that can be studied. Interesting applications should 
be also given beyond the close approximation in collision between black holes ^7,28 
The other point is the meaning of the rescaling in time. We know that a small 
perturbation theory is just a series in G, being G the Newton constant. A dual 
expansion has meaning as an expansion of 1/G", being n a positive integer as we 
will see below. In turn, this means that the scaling of time as t\fG^ in this limit 
implies a prevalence of time variable with respect to other ones into the evolution 
equations. This should be compared with the formal parameter A introduced above. 

It is interesting to note that what we have here is a gradient expansion that 
now appears the dual counterpart to the weak expansion for a gravitational field 
giving us an interesting interpretation of this expansion as the strong coupling series 
for general relativity. 



5. Applications 

In order to show in a given example that a gradient expansion is indeed a 1/G 
perturbation series we consider a perturbation on a Schwarzschild black hole and 
analyze it very near the horizon of the singularity. 

First of all, when one puts [vg = 2GM is the Schwarzschild radius) 

a = ao = 1 - — (33) 
r 

^(0) = 

Irr 



1 - ^ 

r 

(0) 2 (0) 2 • 2 , 

Yac = -r 7Ii = -r sm ( 



into eqs.(32), it is easily verified that higher order corrections are all zero as it should 
be as the Schwarzschild solution is exact. This is a very easy consistency check. 

Now, we change the physical situation by imposing on the lapse function the 
choice 

a{r,t) = 1 - ^ + Ao— sm{ujt) (34) 



that is we perturb the black hole with a radial perturbation of pulsation u) and 
amplitude Aq. This is a standard approach to perturbation in general relativity. 
What we want to check is the behavior with respect to G of the first order correction. 
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The solution of eqs.(32) is straightforward to obtain giving the following solu- 
tions 



1 



1_ is 



,9 1 3 — 4 COS + cos 2wi , / ra\ Lot — smut 



rro 



lee 



1 - 2A 



V r J Tn 



r 3 — 4 cos ujt + COS 2ujt 



+ .. 



2A 



o(l-^) 
V r J r, 



74>4> = lee sin^ 0. 



(35) 



Then, when we approach the horizon of the singularity, that is for r : 
to get 



rg, it is easy 



1 



Irr ~ 2 

1 

lee « -''^ 
= lee sin^ 



1 



^0 



3 — 4 cos ujt + cos 2u>t 



4a; 2 



+ ... 



(36) 



where we recognize that the first order correction goes like l/cj^r^ that is, it is a 
l/G^ term. This is consistent with our expectations given in the preceding section. 

As a final consideration, we know that the gradient expansion at the leading or- 
der admits the cosmological solutions due to Belinski-Khalatnikov-Lifshitz (BKL) 
24,25,26 view of our analysis, these solutions are a generic effect at the leading 
order of a 1/G expansion, and so they appear as a rather general behavior in such 
physical situations supporting recent numerical findings describing the behavior 
of solutions of the Einstein equations near a singularity that seems to be homoge- 
neous in space. 



6. Conclusions 

We have built an expansion for Einstein equations in vacuum for the description of a 
strong gravitational field. The main conclusion is that the spacetime is homogeneous 
near a singularity. An open question is left: What is the proper Bianchi solution in 
the neighborhood of a singularity? The kind of Bianchi spacetime to be considered 
relies in the end on the proper setting of the initial conditions. Two other important 
conclusions to be drawn is that spacetime near the singularity of the big bang was 
homogeneous. A future theory of quantum gravity will have the duty to explain 
such a behavior near a singular point of the classical theory. 

Further studies are deserved when matter is present. In this case the nature of 
the singularity may change "^^ . 

So, we have reached relevant conclusions by the analysis of the Einstein equations 
in vacuum in a strong coupling regime. A lot of physics could be extracted by further 
analysis of this perturbation equations in different situations that at present can be 
managed only numerically. 



1, 2008 21:33 WSPC/INSTRUCTION FILE LANL 



Strong coupling expansion for general relativity 13 

Acknowledgments 

This paper has been inspired by the beautiful lectures on general relativity of 
Vladimir Belinski that I listened at University of Rome on 1992 where I learned 
for the first time about BKL solution and gravitational solitons. 

References 

1. D. Pnueli and C. Gutfinger, Fluid Mechanics, (Cambridge University Press, Cam- 
bridge, 1992). 

2. A. H. Nayfeh, Introduction to Perturbation Techniques, (Wiley, New York, 1981). 

3. C. M. Bender, F. Cooper, G. S. Guralnik and D. H. Sharp, Phys. Rev. D 19, 1865 
(1979). 

4. C. M. Bender, F. Cooper, G. S. Guralnik, H. Moreno, R. Roskies and D. H. Sharp, 
Phys. Rev. Lett. 45, 501 (1980). 

5. C. M. Bonder, F. Cooper, R. Kenway and L. M. Simmons, Jr., Phys. Rev. D 24, 2693 

(1981). 

6. M. Frasca, Phys. Rev. A 58, 3439 (1998). 

7. M. Frasca, Phys. Rev. A 60, 573 (1999). 

8. M. Frasca, Phys. Rev. D 73, 027701 (2006); Erratum-ibid., 049902 (2006). 

9. R. Arnowitt, S. Deser, and C. W. Misncr, in Gravitation: an introduction to current 
research, Louis Witten ed. (Wiley, 1962), chapter 7, pp. 227-265, [gr-qc/0405109]. 

10. M. Alcubierre, gr-qc/0412019. 

11. B. K. Berger, "Numerical Approaches to Spacetime Singularities", Living Rev. Rela- 
tivity 5, 1 (2002). URL (cited on 26 august 2005): http://www.livingreviews.org/lrr- 
2002-1. 

12. H. Friedrich, gr-qc/0508016. An original view of the approach to singularity is given 
in T. Damour, M. Henneaux, and H. Nicolai, Class. Quant. Grav. 20, R145 (2003). 

13. R. Penrose, Phys. Rev. Lett. 14, 57 (1965). 

14. S. W. Hawking, Phys. Rev. Lett. 15, 689 (1965). 

15. S. Hawking, and R. Penrose, Proc. Roy. Soc. A 314, 529 (1970). 

16. D. Garfinkle, Phys. Rev. Lett. 93, 161101 (2004). 

17. R. M. Wald, General Relativity, (The University of Chicago Press, Chicago, 1984), 
pp.183-187. 

18. C. Teitclboim, in Quantum theory of gravity, cd. S. Christenscn (Adam Hilger, Bristol, 
1984) p.327. 

19. R. Jackiw, in Quantum theory of gravity, ed. S. Christensen (Adam Hilger, Bristol, 
1984) p.403. 

20. E. D'Hoker, and R. Jackiw, Phys. Rev. D 26, 3517 (1982). 

21. J. Tcschncr, Class. Quant. Grav. 18, R153 (2001). 

22. H. Lu, S. Mukhcrji, C. N. Pope, and K.-W. Xu, Phys. Rev. D 55, 7926 (1997). 

23. G. B. Cook, and S. A. Teukolsky, Acta Numerica 8, 1 (1999). 

24. I. M. Kalathnikov, and E. M. Lifshitz, Phys. Rev. Lett. 24, 76 (1970). 

25. V. A. Bchnski, I. M. Kalathnikov, and E. M. Lifshitz, Adv. Phys. 19, 525 (1970). 

26. V. A. Belinski, I. M. Kalathnikov, and E. M. Lifshitz, Adv. Phys. 31, 639 (1982). 

27. R. Price and J. PuUin, Phys. Rev. Lett. 72, 3297 (1994). 

28. A. Abrahams and R. Price, Phys. Rev. D 53, 1963 (1996). 

29. D. S. Salopek and J. M. Stewart, Class. Quant. Grav. 9, 1943 (1992). 

30. A. D. Rendall, gr-qc/0503112. 



